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ABSTRACT 


The hydroelastic behavior of very large floating structures (VLFSs) is investigated based on the 
proposed multi-modules beam theory (МВТ). To carry out the analysis, ће VLFS is first divided 
into multiple sub-modules that are connected through their gravity center by a spatial beam with 
specific stiffness. The external force exerted on the sub-modules includes the wave hydrodynamic 
force as well as the beam bending force due to the relative displacements of different sub-modules. 
The wave hydrodynamic force is computed based on three-dimensional incompressible velocity 
potential theory, and the boundary element method with the free surface Green function as the 
integral kernel is adopted to numerically find the solution. The beam bending force is expressed in 
the form of a stiffness matrix. The coupled motion equation is established according to the 
continuous conditions of the displacement and force. The motion response defined at the gravity 
center of the sub-modules is solved by the multi-body hydrodynamic control equations, then both 
the displacement and the structure bending moment of the VLFS are determined from the stiffness 
matrix equations. To account for the moving point mass effects, the proposed method is extended 
to the time domain based on impulse response function (IRF) theory. The accuracy of the 
proposed method is verified by comparison with existing results. Detailed results through the 
displacement and bending moment of the VLFS are provided to show the influence of the number 


of the sub-modules, and the influence of the moving point mass. 
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1. Introduction 


A very large floating structure (VFLS) is a unique type of oceanic structure that embraces a 
range of unprecedented parameters. VLFSs are designed primarily for floating airports and as a 


dock for calm waters on open seas. 

Because it is larger than existing floating structures, A VLFS’s flexibility must be taken into 
account. Also, for a VLFS to be used as an airport or bridge, time-varying moving loads must be 
considered, e.g., the analysis of the displacement and bending moment of a very large floating 
bridge in calm water with some moving point masses (Wu et al., 1998; Fu, 2005), or the 
displacement response of a VLFS in waves with a moving load point (Endo, 2000; Shin et al., 
2001). 

Simulations of the displacement and bending moment of a VLFS are mostly based on 
hydroelastic theory to date. In the frequency domain, hydroelastic theory was initially developed 
in two dimensions (Betts et al., 1977; Bishop and Price, 1977). Wu (1984) extended the work of 
Bishop and Price to three dimensions, and Chen et al. (2003a) developed three-dimensional 
nonlinear hydroelastic theory. Hydroelastic theory can be also extended to the time domain (Wang 
Dayun, 1996; Liu and Sakai, 2002). 

To account for the effects of moving point masses, simulations must be carried out in the time 
domain. Chakrabarti (2001) and Jacobsen (2006) analyzed the hydrodynamic interactions of a 
multi-body system based on the impulse response function method (IRF), but the bodies are not 
connected to each other. Shen et al. (2002) presented procedures for numerical solutions of system 
motion responses of a multi-rigid-body in the time domain. A moving load on a rigid body has 
also been studied (Wu and Sheu, 1996; Shen et al., 2003). 

An approximation theory for the analysis of the displacement and bending moment of a 
VLFS, based on multi-modules beam theory (МВТ), was recently developed by Lu et al.(2015, 
2016). In this method, the VLFS is divided into sub-modules that are connected by a beam with 
specific stiffness through the gravity center of the sub-modules. The external force exerted on the 
sub-modules includes the wave hydrodynamic force as well as the beam bending force due to the 
relative displacements of different sub-modules. Then, the final hydroelastic response equations 
can be established based on multi-body hydrodynamic theory and beam theory. The displacement 


along the VLFS can be obtained by solving these equations, and then the structural deformation as 


well as the section loads can be computed. 

In this paper, Lu's methodology for the computation of the displacement along the VLFS has 
been improved, and some new methods have been developed for the computation of the section 
loads. We discard the beam bending method of Lu et al.(2015, 2016) and calculate the 
displacement response of the structure by the method of matrix transformation. To obtain the 
complete bending moment distribution, the high-order difference is obtained for a few accurate 
moment values. This method is more rigorous, and the displacement and bending moment along 
the whole VLFS can be obtained. This paper avoids the resonance problem caused by small gap 
reasonably when calculating the hydrodynamic coefficient (Miao et al., 2000). 

The methodology has been extended to time domain analysis based on IRF theory, and this 
paper attempts to simulate the hydroelastic response of the VLFS in waves with some moving 
point mass. Study on the moving load of an elastic floating body under wave action is the 
foundation of research on landing pontoons, but there is little literature. This paper combines MBT 
(Lu, 2015), IRF theory (Cummins, 1962; Ogilvie, 1964), and the analysis method of moving load 
on rigid floating body. MBT and IRF are used in the time domain analysis of a VLFS; the floating 
body in the moving load is treated as a rigid body, and the inertia force introduced by the moving 
load is added to the wave force of this floating body. Although the method is very approximate, 


the results can reflect the response characteristics of the structure. 


2. Mathematical formulation 


Lu’s method is based on potential theory and multi-body theory. For the hydrodynamic aspect, 
the ideal fluid assumption is adopted, i.e., the fluid is inviscid, irrotational, and incompressible. 
The incident wave amplitude is assumed to be small relative to a characteristic wavelength and 
body dimension. For the structural aspect, a stiffness matrix is introduced between two adjacent 


modules to consider the flexible structure. 


2.1 Hydrodynamic model 
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Fig. 1 Definition of the fluid and structure boundaries 


Based on the assumptions of an ideal fluid and linearity, the velocity potential can be 


decomposed into three parts as follows: 
P=, + by +x (D 
Where ør, p and фк denote, respectively, the incident wave potential, diffraction wave 


potential, and radiation wave potential.The incident, diffraction potential and radiation potential 


satisfy the following boundary conditions: 
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As shown in Fig. 1, Q is the fluid domain, and Sr, Sg, эк, and 5»аге the free surface, bottom 
surface, wetted body surface of the k" module, and boundary surface at infinity of the fluid, 
respectively. ri, represents the outward-directed unit vector normal to the wetted surface of the Kn 


th 


module. ¢ is the velocity potential. After 


5, 


module. V, is the speed on the wetted surface of the k 


the velocity potential ¢ is obtained, the added mass and the radiation damping of the bodies, as 


well as the wave excitation force, can be calculated. 


4 


In the frequency domain, the excitation forces are related to the incident and diffracted 


potential as follows, 


F., = pio|| (4 + 4.) ri as (3) 
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where J, F 


"а P, 0 and Sog are the imaginary unit ,wave excitation forces, fluid density, wave 
frequency and average wetted surfaces, respectively. 
The added mass and damping coefficient is given by: 
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2.2 Motion equation 
According to Lu’s method, a VLFS is divided into many modules, and the equation of six 
degree-of-freedom (DOF) motion in the frequency domain for a given wave amplitude and wave 


frequency o can be written as follows: 
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where m is the number of floating bodies, [М] is the mass (or inertia moment) matrix of body k, 
[ax] 15 the added mass (or added mass moment) matrix of body k caused by the motion of the 
body itself, [В] is the radiation damping matrix of body К caused by the motion of the body itself, 
[С] is the hydrostatic stiffness matrix of body k, [a,j] is the added mass matrix of body k induced 
by the motion of body j, [b,j] is the radiation damping matrix of body k caused by the motion of 
body j, {fen} is the wave excitation force (or moment) of body k, [U] is the displacement matrix of 
m bodies, [K] is the stiffness matrix that shows the relationship between displacement and force, 
and {их} is the six degree-of-freedom (DOF) displacement of body k expressed as (Xk, Yk, Zk Oks Pro 
У). The number of equations is the same as the number of floating bodies. 


Further, considering all modules, Eq. 5 is extended to the form of an equation set, 
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As a linear system, the coefficient matrix of unknowns [U ] 


„а сап be replaced by [K], and 
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Eq. 6 can be rewritten as 
[E T s Lm m Lone (7) 


where [K7 is the total stiffness matrix. The expressions for the stiffness matrix [K] will be 


discussed in the following section. 


2.3 Analysis of spatial beam 

The six DOF motion of a module's center is restricted by the deformation condition of the 
equivalent beam between two modules. 

As Lu (2016) has written, each side of the beam element i-j has six components of 


displacements and forces: 
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where (Fx; Fy; Fz;) are forces in the i" cross section, and (Мх, My; Mz;) are the corresponding 

bending moments. Similar expressions can be obtained for the j” cross section. Generally, the 

forces and moments are assumed to be imposed on the center of the cross-section of the beam 


element. 


The relationship between [U] and [F] is: 
[F] =[K]., lU], (9) 


where [K];; is the stiffness matrix that is used in Eq. 5, which is shown in the appendix. The 


specific instructions for the matrix are described by McGuire et al. (2000). The size of the matrix 


is 12x12. The matrix is affected by the size of each body, Young’s modulus, and the Poisson ratio. 
When there are m floating bodies, the size of the total matrix [K'] in Eq. 7 will be 6mx6m. 
3. Hydroelastic response of the VLFS 
3.1 Numerical model 

To illustrate this methodology, we use the model mentioned by Yokosuka (Yago and Endo, 
1996). This VLFS is a scaled model of the Mega-Float, constructed and developed for use in 
sheltered waters. The wave amplitude is 1 m, and the main parameters are listed in Table 1. The 
grid in the hydrodynamic calculation of the УГЕ is 160x10x10. 


Table 1. Details of the VLFS model 


Parameters Units Value 
Length L/m 300 
Width B/m 60 
Depth D/m 2 
Draft d/m 0.5 
Young's modulus E/N/m 1.19e10 
Poisson ratio у 0.13 
Density p/kg/m? 256.25 
Water depth H/m 58.5 


3.2 Displacement along the VLFS 

Lu (2015, 2016) proposed a method for the computation of the displacement along the VLFS 
based on spatial beam theory. In this paper, a more generalized method will be developed. 

It can be seen from Eq. 9 that once the displacement [U] has been solved, the section load [F] 
can be computed. The displacement along a given spatial beam can be calculated based on the end 


interfaces' displacements and section loads. 
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The end interface loads F; and F;can be obtained once the displacement И; and U;have been 


Eq. 9 can be rewritten as 


solved. 


To get the displacement U; and section load F; defined at any position between the two end 


interfaces, Eq. 10 is modified as follows: 
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Lu’s method (Lu, 2015) can only calculate the vertical displacement response.Eq.12 can 


(12) 


provide the six-degree-of-freedom displacement along the VLFS, which is a big improvement 


over the original method proposed by Lu (2015, 2016). 

The present results are compared with experimental results (Yago and Endo, 1996) and 
numerically calculated data from three-dimensional hydroelastic theory (Fu et al., 2007), which 
are shown in Fig. 2. In this part, the continuous VLFS is divided into five modules with 0-degree 
wave direction. Values calculated in this paper are denoted by “Present.” Experimental values 
(Yago and Endo 1996) are denoted by “Experiment,” and numerically calculated data from 
three-dimensional hydroelasticity theory (Fu et al., 2007) are denoted by “Fu et al 2007.” These 
results are in good agreement. The hydrodynamic responses of the model in 30- and 60-degree 
wave directions are calculated. The vertical displacement response of the model is compared with 
the results of three-dimensional hydroelasticity theory (Fu, 2005). The results are listed in the 


appendix(Fig.15). 
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Fig.2 Vertical response contrast of VLFS 
3.3 Bending moment distribution along the VLFS 
To obtain the moment of the VLFS, Lu (2016) proposed an equal force method based on the 
Darren Bell principle. It is equivalent to m sets of forces being exerted on the center of each 
module for the VLFS. The whole problem can be described as follows: a free beam is subjected to 


a number of concentrated forces to reach static equilibrium, and to obtain the moment distribution 


of the beam section. This is shown in Fig. 3. The forces conform to the following formulae: 


tae = [К]... [U Jera (13) 
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Fig. 3 Schematic diagram of the equilibrium forces exerted on a free beam 


Lu adopted this method and computed the bending moment along the VLFS, where eight 
sub-modules were used, and the numerical results were compared to the values calculated by the 
hydroelastic theory (Fu et al., 2007). It can be seen from Fig. 4 that the results based on beam 
theory are in good agreement with the results calculated by Fu et al. Since the continuous 
distribution forces are simplified to be the concentrated force on the gravity center of the 


sub-module, only the connection interfaces of the sub-modules are exactly right, which is 


explained in the next section, i.e., eight sub-modules can obtain nine exact points, then the 


bending moment distribution along the VLFS can be computed by B-spline interpolation. 
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Fig. 4 Moment distribution of the model 


To demonstrate the influence of the number of sub-modules, Fig. 5 shows the comparison of 
bending moment distributions of five and eight sub-modules, respectively, in which the frequency 
is chosen as 0.6 rad/s. Here, FBP denotes the six exact values computed by the five sub-modules; 
FPS represents the bending moment distribution along the VLFS via B-spline interpolation; EBP 
means the nine exact values of the bending moment; and EPS represents the bending moment 


distribution along the VLFS via B-spline interpolation. 
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Fig. 5 Moment distributions with different numbers of sub-modules 


It can be seen that there is a big difference between the bending moment distributions for five 
sub-modules and eight sub-modules. To overcome this problem, some improvements have been 
proposed in the following section, specifically to use fewer sub-modules to obtain high accuracy 


of the bending moment distribution. 


3.3 High order interpolation method 

In the multi-sub-module beam theory of Lu (2015, 2016), the continuous distribution 
pressure is treated as the concentrated force exerting on the gravity centers of the sub-modules. 
Therefore, the displacement as well as the bending moment distributions can be determined by 


four parameters, as seen in the following equations: 
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However, the continuous distribution pressure has to be taken into account when solving for 


the displacement and the bending moment, as seen in the following equations: 
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where V(x) is the vertical displacement at x. z,, д, М,, №, M,, and N, are the 


displacement, rotation, equivalent bending moment and shear force, exact bending moment, and 
shear force, respectively, for the center of gravity of the i-th cross section; q(x) is the 
continuous distribution pressure; and x is the distance of the target point from the center of the i-th 
cross section. Eqs. 15 through 18 explain why Lu’s method can only obtain several exact values. 
In Lu’s original work, only the exact bending moment is used to perform the spline 
interpolation. However, it can be seen from eq. 16 that on the exact point, both the bending 
moment value and its partial derivative with respect to the variable x are known. Therefore, the 
higher-order interpolation scheme can be applied. Li (1986) showed that the interpolation curves 


can be cubic order if both the point value and its partial derivative at the two ending points are 
known, i.e., lis Yo yo} and fx, Yo У! } are known. Then, the interpolation curve can be 
expressed as 

p(x)=atbxt+ex + dx (19) 


where 


a jose Ae ОА 
b _ Fox de q (20) 
c D 1 2x% X5 
d 0 1 2x З || и 


Since the current method takes into account both the bending moment and the shear force 
when solving the bending moment distribution along the VLFS, the interpolation accuracy must 
be greater than that of the original method of Lu. Fig. 6 gives the values of the bending moment 
computed by five sub-modules and eight sub-modules, respectively. Here, FBHOI means the 
values are computed by the five sub-modules through the high-order interpolation method, and 
EBHOI means the values are computed by the eight sub-modules through the high-order 
interpolation method. It can be seen that the trend of the values computed by FBHOI is the same 
as that computed by EBHOI, which demonstrates that the HOI method is more accurate than the 
spline interpolation scheme. Besides, the comparison of EBS and EBHOI shows that if the number 
of sub-modules is enough, both the HOI method and the spline interpolation scheme can give 


reasonable results for the bending moment distribution along the VLFS. 
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Fig. 6 Moment distributions computed by the high order interpolation method 


Though the trend of the values computed by the HOI method is the same as for the spline 
interpolation scheme, there is still some difference in the absolute values of the bending moment 
distribution. Therefore the displacement-to-displacement to moment method (DDM) is developed 
to improve the accuracy of the absolute values. Fig. 7 gives the displacement distributions along 
the VLFS computed by five sub-modules and eight sub-modules, respectively; here, the wave 
frequencies are chosen as 0.3 rad/s and 0.6 rad/s. It can be seen that the displacement distributions 
along the VLFS computed by five sub-modules are nearly the same as those computed by eight 
sub-modules. The core idea of DDM is that the displacement distributions calculated by five 
sub-modules and eight sub-modules are considered the same, so the six DOF of eight modules can 
be converted by the six DOF of five modules, and then the bending moment can be calculated by 
eight sub-modules. In the HOI scheme, the bending moment and the shear force are taken into 
account, but the DDM scheme also considers the displacement information, which means this 


method is more accurate than the HOI scheme in application. 
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Fig. 8 Bending moment distribution with different interpolation schemes 


Fig. 8 gives the bending moment distributions along the VLFS computed by the different 
interpolation schemes; the frequencies are chosen as 0.3 rad/s, 0.4 rad/s, 0.6 rad/s, and 0.7 rad/s. 
Here, FTEP means the eight bending moment points computed by the DDM scheme using five 
sub-modules, and FTEHOI represents the bending moment distributions along the VLFS 
computed by the DDM scheme using five sub-modules. It can be seen that the DDM scheme is 


better than the HOI method, and in realistic computations, the combination of five sub-modules 
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and the DDM scheme can produce a reasonable result. 
In the actual situation, considering the computational efficiency, fewer modules are used. The 


DDM method is useful to obtain high-precision calculation results. 
4. Hydroelastic response with moving mass points 


4.1 Time domain hydroelastic theory 

To take into account the effects of moving mass points, the simulation must be carried out in 
the time domain. In terms of the convolution of body motions according to the Cummins 
equations (Cummins, 1962) and Ogilive relations (Ogilvie, 1964), for a multi-body system (with 


N bodies), the equations of motion can be rewritten in the following form: 
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where M and A are the mass and the added mass with infinite coefficient, respectively; B is the 
delay function matrix; C is the hydrostatic resilience matrix; u and the over dot means the 
displacement and the partial derivative with respect to time; and F is the wave exciting force. 

The delay function and the added mass with infinite coefficient can be computed via the 


following equations: 


B,(t)= Sp (@)cos(at)da (22) 
c 0 
M =a(o,,)+ — [na sin(a,,t)dt (23) 
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where a and b are the frequency domain added mass and damping coefficient, and @,, is an 


arbitrarily chosen frequency. 
Based on Lu's theory, for the time domain hydroelastic simulations, Eq. 21 should include 


the stiffness matrix as follows: 
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-dT (24) 


where K is the stiffness matrix coefficient. 
If there exists some other force, e.g., moving point mass, then the above equations should be 


modified as follows: 
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where fis the additional external force. 

Section 3 has shown that the combination of five sub-modules and the DDM scheme can 
provide reasonable results for both the displacement and bending moment distributions along the 
VLFS. Here, indirect time domain simulations are carried out, and the time domain results are 
transformed into the frequency domain. The comparisons of time domain (TD) values with the 
frequency domain (FD) results are shown in Fig. 9. The results are in good agreement, but there 
are some numerical errors. This paper uses the Runge-Kutta method in time-domain iteration. 
Normally, numerical methods have computational errors. In the process of using the IRF method, 
we need to truncate the delay function. This is another numerical error. For a number of floating 


bodies, the truncation term will increase, so the error will increase. The vertical response along the 


VLFS is calculated by Eq. 12, and the [U ] „а erm exists calculation error. When multiplied by a 


6n 


coefficient matrix, the whole displacement response exists calculation error. 
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Fig. 9 Comparisons of the time domain values with frequency domain results 


4.2 Consideration of the moving point mass 


M 


Rl 


Fig. 10 Demonstration of the absolute acceleration of the moving point mass 


To describe the position of the moving point mass, two coordinate systems are introduced in 


Fig. 10. These are the space fixed coordinate system OXYZ and the body fixed coordinate system 


oxyz. With the origin o located at р, =(x,,y,,z,) in coordinate system OXYZ, the position of 


the moving point mass is р, =(x,y,z) in the body fixed coordinate system. Therefore, the 
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position of the moving point mass in the space fixed coordinate system can be written as 


B, B. B, (26) 


If the relative velocity of the moving point mass with respect to the VLFS is taken as V, , and 


keeping connected with the VLFS, then the acceleration of the moving point mass in the inertial 


coordinate system can be computed by the following equation (Shen, 2003; Shi min, 2001): 


а, =а, +а, +2(@xv,) (27) 
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where @ is the rotational velocity of oxyz, а, is the absolute acceleration, а, is the transport 


acceleration, d, is the relative acceleration, and 2(@ху,) is the Coriolis acceleration. 


For the linear case, Eq. 27 can be simplified as follows (Guo-ping, 1989): 
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Therefore, the external force exerting on the VLFS coming from the moving point mass can 


be computed by 
f=-M[T] [T,]ä, -2M [r] (xv, )+[T] Е (30) 


where 
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à, = (3,3, 2,.4,. Fa) (32) 
4.3 Hydroelastic response of the VLFS with a moving point mass 
Here, we use a non- rigorous method to study the hydroelastic response of the VLFS with a 
moving mass point. 
Through the combination of time domain hydroelastic theory (Eq. 25) and the external force 
computation equation of the moving point mass (Eq. 33), the hydroelastic response of the VLFS 


with a moving mass point is investigated. If the moving point mass is located at the i-th 
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sub-module, then the motion equation of the VLFS can be computed by as follows: 


M,, + A, (о) Ев Aw (©) ii, (t) 
Ам (œ) ea " T 

‚| By (t-7) 
F : 

°| By, (t-17) (t) 

C, + К, 0 
+ +4 f(t) 

Cy + К 0 


where f(t) is computed by Eq. (30). 


V 
— 


Module 1 Module 2 Module 3 


Fig. 11 Division of the total VLFS 


Module 4 Module 5 


(33) 


The main purpose of this paper is to discuss the applications of the proposed method for the 


computation of the hydroelastic response of the VLFS, with or without the moving point mass. 


Therefore, only one wave condition (wave length 4 = 420 m) is taken as an example. In the 


numerical simulations, the total VLFS is divided into five sub-modules, as shown in Fig. 11; the 


moving point mass is taken as 10° kg; the velocity of the moving point mass is 5 m/s; and the 


moving point mass begins to go on the VLFS at time 60 s. 


Fig. 12 gives the time history of the displacement and bending moment of the VLFS, while 


Figs. 13 and 14 give the contour of the displacement and bending moment of the VLFS with a 


moving point mass, respectively. It can be seen that before the moving point mass moves on the 


VLES, the time variation of the VLFS is stable, and then the figure shows the influence of the 


moving point mass on the bending moment and the displacement. 
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Fig. 14 Contour of the bending moment of the VLFS 


Conclusion 


This paper first describes a new method, based on both multi-body hydrodynamics and the 
Euler-Bernoulli assumption, to calculate the response of the VLFS. Then two methods, the HOI 
method and DDM method, are proposed to improve the accuracy of the computation of the 
bending moment distributions along the VLFS. Via Fourier transformations, the method has been 
extended to time domain analysis, and then used to investigate the hydroelastic response of the 
VLFS with a moving point mass on it. 

Some test cases have been carried out. All the results are compared with the values obtained 
by three-dimensional hydroelastic theory and the model experiment (Yago and Endo, 1996). The 
results show that the new method is accurate enough to deal with the hydroelastic response of the 
VLES. Some initial attempts have been made to study the hydroelastic response of the VLFS with 
a moving load. Further investigation will be conducted to study the influence of the moving point 


mass on the displacement, as well as the bending moment distributions of the VLFS. 
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Fig. 15. The hydrodynamic response of VLFS in different wave direction 
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